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Abstract 



We study mechanical problems with multiple solutions and introduce a thermodynamic framework 
to formulate two different selection criteria in terms of macroscopic energy productions and fluxes. 
Simple examples for lattice motion are then studied to compare the implications for both resting and 
moving inhomogeneities. 
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1 Introduction 

Many equations describing mechanical systems exhibit more than one solution, and a common problem 
is then to select physically reasonable solutions. The central aim of the paper is to analyse selection 
conditions for problems with inhomogeneities, such as sources or phase interfaces. 

In the literature one finds three commonly cited selection criteria. Two of them are reminiscent 
of Sommerfeld, who gives two physical reasonings of his radiation condition for standing sources 
[Som49, Som62]. He first states that sources have to be sources, and then stipulates that energy 
cannot flow in from infinity. The third selection criterion is the causality principle introduced by 
Slepyan [SleOl, Appendix A], see also the monograph [Sle02], which determines the admissible group 
speeds on both sides of the inhomogencity. All these criteria have been applied mainly in the context 
of Fourier analysis. 

In this paper, we present a different approach to selection criteria, based on thermodynamic fields, 
which are derived via hyperbolic scaling limits and Young measures and describe the underlying me- 
chanical systems on large spatial and temporal scales. In particular, we show that Sommcrfcld's two 
formulations of his radiation condition are naturally related to rather different fields such as produc- 
tion of oscillatory energy and radiation flux. Our framework covers arbitrary nonlinear Hamiltonian 
lattices and PDEs, but we have to specialise to linear or almost linear cases in order to obtain explicit 
expressions for the thermodynamic fields. 

While we recover Sommcrfeld's results for standing sources, the case of a moving phase interface 
turns out to be remarkably different. Specifically, in the first case we find that both thermodynamic 
selection criteria are equivalent and agree with the causality principle. In the second case, however, the 
condition on radiation (energy) fluxes can be in contradiction to the entropy inequality, the production 
criterion and the causality principle. Thus the energy flux criterion should not be applied uncritically 
to problems with moving inhomogeneities. 
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Wc present our exposition along two guiding examples of nearest neighbour (NN) chains of atoms. 
The governing equations are 

Xj(t) = &(x j+1 (t) - Xj (t)) - $'(*,•(*) - Xj -i(t)) + at)5 j0 , (1) 

where £ is the external forcing of atom j = 0, and the two examples are (i) a harmonic chain with 
periodic forcing £, (ii) a moving phase interface in a chain with bi-quadratic interaction potential $ 
and no external forcing. 

Sommerfeld's radiation condition Sommerfeld's approach to radiation conditions is described in 
his book [Som62, §28], see also [Som49], and we thus just describe the gist of some key arguments. 
However, the discussion given in §3 for the case of a lattice model closely resembles the original spatially 
continuous case, and we refer the reader to that section for a more detailed mathematical treatment. 

Sommcrfcld considers the wave equation on M 3 with a temporally periodic forcing at the origin. 
A separation of variables ansatz leads to the Hclmholtz equation with a source at the origin, thus 
an inhomogeneous equation. The corresponding homogeneous equation has nontrivial bounded ker- 
nel functions (in M 3 ). Some kernel functions are excluded by boundary or symmetry conditions, in 
Sommerfeld's case by requiring that the solution is radially symmetric and decaying at infinity. This 
still does not single out a unique solution. In particular, there are two singular and decaying solutions 
u^ 1 to the inhomogeneous Hclmholtz equation, which are well-defined outside the origin. Let x^ be 
the corresponding solutions of the forced wave equation. The choice of the direction of time made in 
the separation of variables ansatz renders the corresponding solution x + an outwardly radiating radial 
wave and x~ an inwardly radiating wave. In fact, x + describes a source solution to the forced wave 
equation in the sense that the forcing supplies energy to the system. On the contrary, x~ corresponds 
to a sink solution, as the forcing now deprives the system of energy. 

Sommcrfcld now introduces a binary choice, allowing only waves which propagate outwards and 
dismissing those which propagate inwards. This selection is necessary both for physical and mathe- 
matical reasons: mathematically, the choice is an integral part of the arguments leading to a unique 
fundamental solution of the forced wave equation. Physically the two solutions x^ are qualitatively 
very different. In Sommerfeld's words [Som62], 

"Qucllcn sollen Quellen, nicht Scnkcn dcr Energic scin." [Sources have to be . 
sources, not sinks of the energy] 

We call this the first formulation of Sommerfeld's radiation condition. He then gives what we call the 
second formulation, 

"Die von den Quellen ausgestrahltc Energie mufi sich ins Uncndliche zerstreucn, 

Energie darf nicht aus dem Unendlichen in die vorgeschriebenen Singularitaten 

des Feldes eingestrahlt werden. v [The energy radiated from the sources has to (SOM2) 

scatter to infinity, energy must not be radiating from infinity into the prescribed 

singularities of the field.] 

Thermodynamic approach In §2 we establish a thermodynamic framework to analyse (SOM1) 
and (SOM2). The starting point are discrete local conservation laws, which can be derived from (1) 
and describe the transport of mass, momentum and energy on the microscopic scale. In a first step 
we introduce a small scaling parameter e and use the concept of Young measures to establish the 
macroscopic limit e —> 0. In this limit the microscopic conservation laws converge to their macroscopic 
counterparts, which arc PDEs with densities, fluxes and productions. 

In a second step wc identify the macroscopic quantities that admit a thermodynamic interpretation 
of Sommerfeld's radiation condition. The radiation flux Q is defined as the Galilean invariant part of 
the energy flux and (SOM1) is naturally related to the sign of Q on both sides of the inhomogencity. 
We also introduce the oscillatory energy E osc , which measures the amount of macroscopic energy 
stored in microscopic oscillations, and the non-oscillatory energy E non . Sommerfeld's first formulation 
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of the radiation condition can now be reinterpreted in terms of energy productions. More precisely, 
for forced excitation problems as studied by Sommerfcld for the wave equation, (SOM1) is naturally 
related to the sign of the local production of total energy E = E osc + E non due to the external forcing. 
The situation is different for phase transition waves as these are not driven by external forces and 
conserve the total energy exactly (in a local sense). Instead, the key energetic phenomenon is now 
that the propagating phase interface causes a constant transfer between oscillatory and non-oscillatory 
energy. This transfer gives rise to a production of oscillatory energy, and (SOM1) is naturally related 
to the sign of this production. 

Notice that the above arguments are solely based on energy productions and fluxes, that means on 
fundamental concepts from continuum mechanics and thermodynamics. They can therefore, at least in 
principle, be also applied to more general systems, as for instance NN chains with generic double-well 
potentials. However, only in some special cases it is also possible to derive explicit expressions for all 
densities, fluxes and production terms. 

We finally emphasise that our macroscopic interpretations of (SOM1) and (SOM2) provide only 
necessary conditions. To ensure uniqueness, in particular to exclude linear combinations of x + and 
x~ in Sommerfeld's forced excitation problem, one typically imposes microscopic selection criteria, 
which also exploit non-macroscopic properties of solutions. Examples are the asymptotic radiation 
condition (21) or, equivalently, the causality principle. Both methods are inextricably connected 
with Fourier methods. It is therefore not clear how to generalise these selection criteria to nonlinear 
problems. 

Results for the special examples In §3 we apply the thermodynamic framework to the analogue 
of Sommerfeld's problem in the harmonic NN chain and find that (SOM1) and (SOM2) arc equivalent 
and comply with the microscopic selection criterion. Afterwards we study phase transition waves in 
§4 and show that both conditions provide different selection criteria for moving phase interfaces. We 
show that (SOM1) is equivalent to the usual entropy inequality, and thus we propose to apply (SOM1) 
to problems with moving inhomogeneities as well. The application of (SOM2), however, turns out 
to be problematic. Specifically, we show that for subsonic wave speed close to the speed of sound, 
the selection criterion (SOM2) selects a phase transition wave that violates the entropy inequality as 
well as the causality principle. We thus conclude that (SOM2) has to be discarded, at least for some 
problems with moving inhomogeneities. 

2 Macroscopic field equations in the presence of microscopic 
oscillations 

Our approach to Sommerfeld's radiation condition in §3 and §4 is based on macroscopic balance laws 
that govern the effective dynamics of averaged quantities on large spatial and temporal scales. In this 
section we derive and discuss these balance laws and provide the atomistic expressions for all densities, 
fluxes and production terms. Our exposition is formal but we emphasise that all arguments can be 
made rigorous using Young measures; sec [DH08, DHR06] and Appendix A. For the sake of clarity, 
we start with the conservation laws for the unforced NN chain, that is (1) with ( = 0. 

Microscopic conservation laws We can rewrite (1) with £ = in terms of the atomic distances 
(discrete strains) rj := and velocities Vj := ij as first order equations 

fj = v j+1 - vj , v 3 = $' fa ) - #'fa_i) , (2) 

which can be viewed as the discrete counterparts of the local conservation laws for mass and momentum 
in Lagrangian coordinates, see (6). We note that (2) constitute local conservation laws in the sense that 
the time derivate of some atomistic observable equals the discrete divergence of another observable. As 
a consequence we find, assuming appropriate decay conditions for j — > ±oo, that each local conservation 
law implies a global one by summing over j. The discrete divergence structure is also an important 
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ingredient for the thermodynamic limit because it allows us to derive the macroscopic conservation 
laws, sec Appendix A. 

Since the unforced NN chain is an autonomous Hamiltonian system with shift symmetry, we can 
also derive a local conservation law for the energy, namely 

f t {\v 2 + Sfa-O) = v^'{r 3 ) „,•_!#'(»•,_!). (3) 

To characterise the thermodynamic properties of NN chains, we now derive a macroscopic description 
by applying the hyperbolic scaling of space and time. For a given scaling parameter < e <C 1, we 
define the macroscopic time t and the macroscopic particle index £ by 

r = et, f = sj, (4) 

but we do not scale distances and velocities. We then regard the atomistic data that correspond to a 
solution of (2) as functions r £ and v e that depend continuously on r and are piecewise constant in £. 
That is, we identify rj(t) = r £ (et,ej) and Vj(t) = v £ (et,ej). 

Of course, the functions r £ and v e will, in general, be highly oscillatory with wave length of order e 
and do not converge as e — > in a pointwisc sense. However, as long as the solution to (2) is bounded 
we can assume, thanks to weak compactness, that for any atomic observable ip = ip{r,v) the functions 
tp{r e ,v e ) converge weakly as e — > 0, see Appendix A. The limit function (ip) is then non-oscillatory 
and can be regarded as the thermodynamic field of ip, which means that (tp}(T, £) gives the local mean 
value of -0 in the macroscopic point (r, £) . 

Macroscopic conservation laws In the thermodynamic limit e — > the discrete conservation 
laws (2) and (3) transform into 

d T {r)=ds(v), d T (v)=d e (&(r)), d T (±v 2 + <P(r))=ds(v&(r)). (5) 

These PDEs describe the local conservation laws for mass, momentum and energy on the macroscopic 
scale and are well known within the thermodynamic theory of elastic bodies. In fact, they can be 
written as 

d T R-d e V = Q, d T V + d s P = 0, d T E + d^F = 0, (6) 

with macroscopic strain R = (r), macroscopic velocity V — (v), pressure P = -($'(r)>, total energy 
density E = {\v 2 + ^(r)), and total energy flux F = —(v$f'(r)\. Moreover, splitting off the Galilean 
invariant part from both the energy density and the energy flux, we find 

E=\V 2 + U 1 F = VP + Q, (7) 

with internal energy density U and heat flux Q. Radiation in the sense of Sommcrfcld precisely means 
energy transport via Q, see also §3, and therefore we call Q the radiation flux. 

We emphasise that, in general, the conservation laws (6) do not constitute a closed system, but must 
be accompanied by closure relations. Unfortunately, very little is known about the thermodynamic 
limit for most initial data and general interaction potentials. In some cases, however, it is possible to 
solve the closure problem. Below we show that all thermodynamic fields can be computed explicitly for 
(i) Sommerfeld's fundamental solution in forced harmonic NN chains and (ii) phase transition waves 
in NN chains with bi-quadratic potential. 

Oscillatory energy For our purposes, the split (7) is not sufficient; it is convenient to introduce 
another split of the total energy density E that accounts for the fact that the computations of local 
mean values and nonlinearitics do not commute in the presence of oscillations. To obtain a precise 
measure for the strength of the oscillations we write 

E = E non + E osc , E non = iy 2 + $(-R), 
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which means 

E non = \(v) 2 + $«r», E osc = §<(«- («)) 2 ) + (*(r) - S«r») (8) 

and furthermore implies {/ = 3>((r)) +£' osc . We refer to -E sc and -E n on as oscillatory and non-oscillatory 
energy density, respectively, and emphasise that E osc measures precisely the amount of macroscopic 
energy that is locally stored within the oscillations. From (5) and (8) we now conclude that the partial 
energies are balanced by 

8 T E osc + dsQ = S, d T E non + 8 ( (PV) = -S, (9) 
where the production terms are given by 

E = -(P + &(R))dtV = (<*'(r)) - &((r)))8s(v) . 

The physical meaning of S becomes apparent when we compute the temporal change of the oscillatory 
and non-oscillatory energy contained in a control volume Q = [a, b\. This gives 

A J" E osc d£ + Q\fl b a = jf 3d£ = - A jf E non d£ (PV)\IZ b a . 

Since the flux terms Q and PV are related to the exchange of energy with the exterior of f2, we conclude 
that E(t, £) measures precisely how much non-oscillatory energy is transferred into oscillatory energy 
at time r at the point £. 

We further mention that the second law of thermodynamics is, at least on an intuitive level, related 
to a sign condition for 3. While each scaling limit of the unforced NN chain conserves the total energy 
according to (5)3, there exist many solutions that yield a non-vanishing field S in the limit e — > 0. 
For S > these solutions 'dissipate' non-oscillatory energy (at the expense of increasing oscilllatory 
energy); this is in perfect agreement with our physical intuition since oscillations can be regarded as 
microscopic fluctuations. If time is reversed, however, 5 becomes negative and the solution extracts 
energy from the fluctuations, which contradicts the intuition. This is discussed in more detail in §3 
and §4. 

Thermodynamic fields for harmonic oscillations With the central quantities E osc and E non 
now being defined, we can return to the analyis of thermodynamic fields. Motivated by the discussion 
in §3 and §4, we restrict the analysis and compute the thermodynamic fields for travelling waves in 
harmonic NN chains. A travelling wave for the NN chain is an exact solution to (1) with £ = that 
satisfies 

Tj (t) = TZ(j - c ph i), Vj (t) = V(j - c ph t) (10) 

for some phase speed c p h and profile functions 1Z and V that depend on the phase variable <p = j — c p hi- 
Travelling waves in NN chains arc determined by advance-delay differential equations, see [FV99, 
DHM06, HcrlO], and describe fundamental oscillatory patterns. 

The hyperbolic scaling (4) implies that each travelling wave solution for the harmonic NN chain 
(10) converges as e — > to a Young measure that is constant in space and time. The thermodynamic 
fields of all observables ip are therefore independent of (r, £ ) and can be computed from R and V as 
follows 

(^) := i im _L f ^(^VfeOJdp. (11) 

L->oo 2L J_ L 

For a harmonic NN chain, which has the interaction potential $(r) = \c\r 2 + d\r + do, we im- 
mediately verify by Fourier transform that for prescribed c p h with < |c p h| < cq travelling waves are 
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given by 



M M 



K(tp) = R + Aj cos (nj(p + m/2 + rg), V(tp) = V =F Co y~] Aj cos (w^y + ffo), (12) 

i=l i=l 

with "— " and "+" for left and right moving waves, respectively, that is, for c p h < and c p h > 0, 
respectively. Here the wave numbers Ki, i = 1, . . . , M, denote the positive solutions to Cp h fc 2 = £l(k) 2 , 
where Q(k) = 2cosin(fc/2) is the dispersion relation of the harmonic NN chain. In particular, near 
sonic waves with c p h ~ cq have M = 1 and depend, up to the phase shift rji, on the four independent 
parameters R, V, A — Ai 7 and k = K\. 

The thermodynamic fields for such harmonic travelling waves can easily be computed by (11) and 
(12). In fact, thanks to (r) = R and (v) = V, wc find 



P = -c£R-d 1 , E non = ±V* + ±4R' + d l R + d , (13) 



as well as 



A. 



E osc = ±c 2 A 2 , Q = ±(c ^^cos(K i /2))s osc (14) 



with A 2 = -^f ■ For periodic waves with M — 1 we therefore have Q = c gI E osc . where c gr = 

± |f2'(«;)| is the group speed. 

We emphasise that the thermodynamic computations presented above can be extended to super- 
positions of finitely many harmonic travelling waves, with obvious modifications. We also mention 
that a complete characterisation of the energy transport in harmonic lattices can be derived in terms 
of Wigner-Husimi measures [Mie06] . 

Macroscopic description of forcing In order to generalise the formalism from above to forced 
NN chains we assume, for simplicity, that the forcing acts only in the particle j = 0, see (1), and that 
C is periodic with 

C(*) = C(* + tper), / P " C(t) d* = 0. (15) 

Jo 

These conditions guarantee that the forcing does not contribute to the macroscopic conservation laws 
for mass, (2)i, and momentum, (2)2. The forcing, however, in general supplies some energy to the 
system, and hence the conservation law (2)3 must be replaced by 

d T E + d 6 F = 6{T)5o(d£). (16) 

The macroscopic energy production at £ = can be computed either as the jump of the macroscopic 
energy flux at £ = or by averaging the microscopic energy production. This reads 



(r) = F(r, 0+) - F(t, 0-) = lim lim - / v (t)((t) dt. (17) 

2d J(r-S)/e 



3 The forced harmonic NN chain 



Here we present the analogue to Sommerfeld's classical problem in harmonic NN chains, that is, the 
localised forced excitation problem 

(d!-c 2 A 1 )x j (t) = at)6° i , (18) 
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where Ai denotes the discrete Laplacian Ai Xj . — Xj-\.\ I X j — \ — 2xj . For simplicity we normalise the 
speed of sound to c.q = 1 and assume that the chain is periodically forced at one of its eigcnfrcqucncies 
a with < er < 2. The analysis presented in this section resembles that of the spatially continuous case 
studied by Sommerfeld [Som62] ; the purpose of this section is to familiarise readers with our approach 
to selection criteria, and to show that we recover Sommerfeld's classical results for a standing source 
in a spatially discrete medium. 

Explicit solutions via Helmholtz equation The separation of variables ansatz Xj(t) = 
Re(uj exp(— \<rt)) transforms (18) into the discrete Helmholtz equation 

a 2 u ] + AiUj = 5 , (19) 

which can be solved by Fourier transform. There exist two special solutions u + and u~ defined by 

± _ L exp(±i/t \j\) 

J i2fi(/s)n'(«) ' 

where re = re(cr) denotes the unique solution to 

cr 2 = f2(fc) 2 , n(jfe) =2sin(fc/2), < k < tt. 

Of course, the special solutions u~ and u + can be linearly combined and also superimposed by 
plane waves with wave numbers ire, which are the kernel functions of the discrete Helmholtz operator. 
The general solution to (19) can therefore be parameterised by a,/3 £ C as 

u.j = Uj(a,(3) = u~j +aexp(— inj) + (3 cxp (+inj). (20) 

Note that in particular = Uj(0, 0) and uj = Uj(a~, /3~) with oT = j3~ = — (i20(re)f2'(re)) . 

Sommerfeld's approach to the radiation condition can be viewed as the endeavour to remove the 
non-uniqueness and to single out a unique choice for a and f). To this end, he introduces a microscopic 
selection criterion, whose analogue in the harmonic NN chains reads 

lim i^p- - inuj) = 0, lim ( ^ + iKUj) = 0, (21) 

j-*+oo \ dj J i^-oo V dj J 

and implies that a = (3 = in (20). In particular, the microscopic radiation condition selects u + but 
rules out u~ . 

Macroscopic aspects of Sommerfeld's radiation condition We now show that the binary choice 
between u~ and u + can be understood in terms of purely macroscopic conditions on the production 
of oscillatory energy and the direction of the radiation fluxes. As remarked in §2, the thermodynamic 
framework can be extended to superpositions of harmonic waves. It is thus possible to characterise all 
bounded solutions to (18), in particular kernel functions and their superpositions. The result of such 
an analysis is that the thermodynamic interpretation of (SOM1) and (SOM2) rejects superpositions 
of waves as long as their influx contribution exceeds the outward contribution. Thus, a half space of 
all bounded solutions is rejected, and a half-space accepted. We show this analysis here in detail for 
the two extreme cases corresponding to «_ and u+. The analysis for the other solutions is, mutatis 
mutandis, analogous yet more complicated terms arise. 

At first we notice that u + and u~ correspond to the real-valued displacements 

± sin(re|j| TQ-t) 

j [ ' 2fi(re)Q'(re) ' 

Each of these solutions to (18) consists of two counter-propagating travelling waves that are glued 
together at j = 0, where the travelling waves propagate towards and away from the inhomogeneity at 
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periodic forcing 




Figure 1: Sommcrfcld's source solution for the harmonic NN chain: The energy pumped in by the 
periodic forcing is radiated towards both +00 and — 00. If time is reversed the source becomes a sink 
and the radiation fluxes on both sides change their sign. 

j = for x + and x~ , respectively. We also notice that x + and x~ transform into each other under 
time reversal, and that they define the atomic distances and velocities 

vf(t) = TAcos{K\j\To-t), rf{t) = Acos{n\j + \\To-t) (22) 

where the amplitude A > is given by 1/A = 2J7'(k) = 2 cos (k/2). 

To determine the thermodynamic limit e — > 0, the key observation is that both x~ and x^ generate 
Young measures that are (i) independent of the macroscopic time r = st, (ii) constant for £ < and 
£ > 0, where £ = ej is the macroscopic particle index, and (Hi) generated by periodic travelling waves. 
These assertions follow directly from (22) and the definition of Young measure convergence, sec §2 
and Appendix A. They also imply that the macroscopic conservation laws (6) are trivially satisfied 
for £ 7^ 0. The relevant macroscopic information is therefore encoded in jump conditions for the 
thermodynamic fields at £ = 0. 

Using (13) and (14) we now conclude that almost all thermodynamic fields are globally constant 
with 

R = 0, V = 0, P = 0, F = Q, E non = 0, E = E osc , E osc = \A 2 , 

while the radiation flux Q is piecewise constant, Q(t, £) = ±sign(£)f2'(/t).E l osc =: Q S gn({)oo- The two 
values for Q are given by 

Q~oo = Tj^, Q + oo = 

and computing the macroscopic energy production by averaging, see (17), we find 

2n/a 2n/a 

0=—£- [ v {t) cos (at) dt = [ cos(at) 2 dt = ±\A. 

2tt J 2?r J 



Sommerfeld's first condition (SOM1) is naturally related to the sign of 8. For x + we have 9 > 0, so the 
forcing pumps in energy at j = and the solution describes a source of the energy. The solution x~, 
however, corresponds to a sink as 8 < implies that energy flows out constantly at j = 0. Moreover, 
for the solutions at hand the balance of total energy reduces to 

8 — Q+oo Q—oo ^(^-i-oq 2Q—OQ, 

and implies that Sommerfeld's first and second formulation of the radiation condition are equivalent. 
Namely, energy that is pumped in at j = must be radiated away, and hence the radiation fluxes 
must point away from j = , see Fig. 1. Conversely, energy that is deprived from the system at £ = 
must be radiated in from ±00. We thus conclude that both (SOM1) and (SOM2) select the solution 
x + but reject x~ . 

We close this section by mentioning that the causality principle, sec [SlcOl, Slc02] and the discussion 
in §4, requires that all oscillatory modes in front and behind the interface satisfy c gr > c p h and 
c gr < Cph, respectively. It is not hard to see that this principle also favours x + and, more generally, 
implies (SOM2) for a standing source. 



8 



Co 
Cl 

C2 

31 271 371 k 

Figure 2: Phase transition waves with |c p h| > C2 have periodic tails and come in two different types. 
Type-I waves have either < c gr < c p h or > c gr > c p h, whereas type-II waves correspond to either 

C gr < < Cph Or Cg r > > c pri . 

4 Phase transition waves for a bi-quadratic NN chain 

We now consider travelling waves with a moving inhomogeneity and apply the selection criteria (SOM1) 
and (SOM2) to these waves. We mention, however, two important caveats of our analysis: firstly, we 
assume the existence of subsonic travelling waves with a single inhomogeneity (interface). Guidance for 
the existence can be taken from [TV05]; yet existence is a subtle issue, and a rigorous existence proof 
is available only for a small regime of subsonic velocities [SZ09, SZ]; it is also worth mentioning that 
there is a velocity regime where no travelling waves with a single interface can exist [SSZ]. Secondly, 
the selection criteria that result from the thermodynamic interpretation of (SOM1) and (SOM2) are 
necessary but not sufficient. 

We now study heteroclinic solutions to the travelling wave equation (10). To calculate the ther- 
modynamic fluxes explicitly, we restrict our considerations to the NN chain with piecewise quadratic 
interaction potential 

$(r) = ±min{(r- l) 2 , (r + 1) 2 }, (23) 

but mention that our thermodynamic arguments can, at least in principle, be generalised to genuinely 
nonlinear potentials as well. (The double- well nature of $ describes the co-existence of different stable 
states and thus the possibility of interfaces between those states.) 

The potential (23) is normalised to have unit sound speed, Co = 1. As illustrated in Figure 2, there 
is a critical velocity C2 > such that for all with c prl with C2 < |c prl | < 1 there is a unique solution 

K > to 

c 2 h fc 2 = n 2 (k), fi(fc) = 2 sin (jfc/2). (24) 

From now on we solely consider waves with |c pri | > C2 because then the tails are periodic with a unique 
wave number k as chosen above. This means, for any travelling wave with a single interface we have 

(11, V)M (K ±00 , V ±00 )(y), (25) 

where both (1Z +001 V+oo) and (TZ-oo, V_oo) are periodic travelling waves (possibly constant) with phase 
speed c pri and group speed c gr . To compute the thermodynamic fields explicitly, it is necessary that 
the asymptotic microscopic strains are confined to the harmonic wells. We thus require that both 
IZ-oo and IZ+oo have a definite sign. By symmetry we can assume that 7£±oo(y) *c 0, and by shift 
invariance we can also assume that TZ(0) = 0. Thus, the interface moves along j = c p ht and £ = c p hT 
in the microscopic and macroscopic space-time coordinates, respectively. Notice, however, that we 
have not fixed the sign of c p h, so the wave travels from negative strain to positive strain for c p h > 0, 
and the other way around for c p h < 0. 

Macroscopic constraints for phase transition waves Under the assumption that travelling 
waves with a single interface as described above exist, all thermodynamic fields are constant on the 
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left and on the right of the interface and are completely determined by the periodic tail oscillations 
in (25). The macroscopic conservation laws therefore reduce to jump conditions via d T ~~> — c p h[] and 
®£ [], so the PDEs (6) transform into 

c P h[P] = -m Cph[Vl = [P], c ph {E aon + E oac ] = lPV + Ql (26) 

Note that the asymptotic jump and mean value of any thermodynamic field X are given by 

(Xj =X +00 -X. 00l {X} = ±(X +00 +X. 00 ), X ±OQ = lim X{t,0- 

We now express the asymptotic values of all thermodynamic fields X in terms of R± cxn V±oo, A± OD , 
and the speeds c p h and c gr . In this way, we recover well-known jump conditions and kinetic relations 
for phase transition waves [Tru82, Tru93]. The strategy of computing thermodynamic quantities as 
averages of atomic observables is well established, see for instance [TV05, SCC05]. However, it is 
usually not based on Young measures and hyperbolic scaling limits. It further seems that the concepts 
of oscillatory and non-oscillatory energy have not been used before in the context of phase transition 
waves. 

Due to (23) and the sign choice for 1Z± 00 , we have P±oo = (R±oo) = —R±oo ± 1 and hence 
[P] = 2 — [i?]. The jump conditions for mass (26)i and momentum (26)2 thus imply 

M = ~r~~ji~^ in - -y^fr, (27) 

and therefore 

Cnh {R} - Cph {V} c ph {R} - C 2 {V} 

[P„o„ J = 2 phl _ 2 P , iPVj = 2 P 1 / 2 phl . 

1 C ph C ph 

Using this and the formulae for P osc and Q from (14), we then find 

-Cphl^non] + [P^] =2Cph{P}, -Cph[P OS c] + [Q] = (C gr -C ph )[P osc ], (28) 

which is the analogue to (9). Consequently, the jump condition for the total energy (26)3 enforces 
that the productions for oscillatory and non-oscillatory energy cancel via 

5 = -2cp h {R} = (c gr - c ph ) [P osc ] . (29) 

This formula is important as it reveals that for phase transition waves there is no production of total 
energy but instead a steady transfer between the oscillatory and the non-oscillatory contributions of 
the energy. This transfer has power —2 {R} and drives the wave. More precisely, the configurational 
force T satisfies 

c ph T = H, T := I$(P)] - {$'(#)} {R\. 

This is the kinetic relation and follows from (27) and (28) thanks to \E non \ = [Vj {V} + [$(#)], 
{PVj = {P} fV] + [P] {V} and P = -$'(P). The production of oscillatory energy S is the process 
commonly called dissipation. Recall, however, that the total energy E = E osc + P non is not dissipated 
but conserved according to the energy laws (6)3 and (26)3. 

We finally notice that time reversal changes the sign of c p h, c gr , S, V±oo, Q±oo, but does not 
affect P±oo, ^4±oo, Pose, ±00, Pnon, ±00, Ptoo, or T. We also observe that all thermodynamic fields are 
completely determined by 

c p h, A-oo, A +QO , {V}. (30) 

In fact, from (30) we compute k and c gr by (24) and set [P OS c] = ^I^ 2 !- Afterwards we solve (27)i 
and (29) for P_oo and P_oo, which then allow us to compute V-oo and V+oo from (27)2- 
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Type-I wave with < c gl < Cph 
radiation flux Q__ "*''\/^*•^"'"\/'' V, • 




Type-II wave with c gl - < < c p h 
radiation flux 




radiation flux Q +a 
wave speed c p h 



Figure 3: Phase transition waves are driven by a constant transfer between the oscillatory and the 
non-oscillatory energy and the radiation fluxes on both sides of the interface have the same sign. The 
cartoons illustrate the source solutions for type-I and type-II waves, which represent different order 
relation for the group velocity c gr and the phase velocity c p h- If time is reversed, the interface becomes 
a sink of oscillatory energy and the radiation fluxes on both sides change their sign. 



Type-I wave according to the causality principle 

Q-„>0 




e+» = o 

Cph > 



Type-II wave according to the causality principle 

2— > o 




e + »,=o 

Cpi, > 



Figure 4: The causality principle selects phase transition waves that propagate into a region without 
oscillations. These waves are, up to Galilean invariance, uniquely determined by their speed c p h- 



The jump conditions derived in this section constitute macroscopic constraints which are necessary 
for the existence of a phase transition wave with speed |c p h| £ (c2, 1). However, it was proven in [SZ] 
that these conditions are also sufficient, at least for near sonic speeds with < 1 — |c p h| « 1. In 
conclusion, there exists a four-parameter family of candidates for phase transition waves, that means 
of hetcroclinic travelling waves. It is now very natural to ask which of them are physically reasonable, 
and so selection criteria come into play. 

The macroscopic aspects of Sommerfeld's radiation conditions Wc first consider (SOM1). 
It is reasonable to require that the interface is a source rather than a sink of oscillatory energy. The 
production S therefore has to be non-negative, which means 

3= (cgr -Cp h )[£ 08c ] > 0. (31) 

This inequality is equivalent to 

CphT > 0, (32) 

which is the usual entropy condition for phase transition waves (see for example [TV05]). For all 
waves considered here, (31) implies -E sc,+oo < -Eosc, -oo for waves moving to the right and E OSCi+oc > 
E sc, -co for left-moving waves. In other words, in all cases we have S > if and only if the oscillations 
have smaller amplitude in front of the interface than behind the interface. (SOM1) select these solutions 
but rejects waves that travel from regions of low oscillations into regions of high oscillations. Note, 
however, that oscillations in front of the interface are not ruled out since it is only required that the 
wave propagates in direction of decreasing oscillations. This implies that there is still a four-parameter 
family of phase transition waves which satisfy (SOM1). 

Sommerfeld's second formulation (SOM2), which stipulates that "energy is carried away from the 
interface" , translates directly into a condition on the radiation flux. It requires, on both sides of the 
interface, that Q points away from the interface. This condition is very restrictive for phase transition 
waves with periodic tails because both Q-oc and Q+co have the same sign as the group velocity c gr , 
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see (14). Thus (S0M2) can only be satisfied if there are no oscillations on one side of the interface. 
However, which side of the interface is chosen by (SOM2) depends on the sign of c gr , and therefore 
we distinguish between two types, see Figures 2 and 3. Type-I waves have c\ < |c p h| < 1, where 
ci := 2£l(n/2)/ir, and this implies sgnc gr = sgnc p h and |c gr | < |c p h|- Type-II waves correspond to 
c-2 < |c p h| < ci, which means sgnc gr ^ sgnc p h- 

For type-I waves, the radiation fluxes behind and in front of the interface point towards and away 
from the interface, respectively. This is illustrated in Figure 3, and holds regardless whether (31) is 
satisfied or not. Energy is therefore always radiated towards the interface, and the second formulation 
of the radiation condition can only be satisfied if there are no oscillations behind the interface. Those 
waves, however, are usually regarded as unphysical as they violate (31) and (32). The only solution 
candidates that would be accepted by both formulations have no oscillations, neither in front nor 
behind the interface; however, such single transition waves do not exist for the potential (23). 

The discussion is different for type-II waves. There is still radiation into the interface but now 
the radiation flux impinges from ahead of the interface. Therefore, both (SOM1) and (SOM2) are 
simultaneously satisfied by type-II waves that propagate into a region without oscillations, that is, 
E OSCt +oo = for right-moving waves and -E OS c.-oo = for left-moving waves. In particular, there 
exists a three-parameter family of type-II waves that satisfy both (SOM1) and (SOM2). 

In summary, for type-I waves, and hence for near sonic waves, (SOM1) and (SOM2) contradict 
each other and would, if applied together, reject any bounded phase transition wave with a single 
interface. For typc-II waves, however, (SOM2) implies (SOM1) and allows for a three-parameter 
family of waves. This is in stark contrast to the case of a standing source discussed in §3, where both 
criteria are equivalent. 



Microscopic selection criteria Besides macroscopic criteria as described above, there also exist 
microscopic selection principles for phase transition waves. These are far more restrictive and select a 
two-parameter family of phase transition waves, as shown below. For the sake of comparison we now 
summarise the main arguments leading to microscopic selection criteria for phase transition waves in 
bi-quadratic NN chains and refer to [TV05, CCS05] for more details. The key idea is that under the 
condition sgaTZ((p) = sgny), each phase transition wave is determined by the affine advance-delay- 
differential equation 

cl h dpl= Aft-Asgn^. 

This equation can be regarded as the analogue to the inhomogeneous Helmholtz equation (19), and 
solutions can be represented by 1Z{tp) = Rc(S((p)) with 

i f il(k) 2 exp(+ikip) 



S ^ = -«j -0**3 dfc ' (33) 

r p 

where V is an appropriately chosen contour in the complex plane. The microscopic selection criterion 
is based on the causality principle [SleOl, Sle02] and requires that all oscillatory modes in front and 
behind the interface satisfy c gr > c p h and c gr < c p h, respectively. The contour V is therefore chosen 
as the dented real axis that passes the origin k = from below but the other real- valued poles of 
the integrand in (33) from above. Jordan's lemma from complex-valued calculus then provides the 
following expressions for the thermodynamic fields for a right moving wave 

i? ±00 = ±-^V + ^^, A . ao = -^- l A +oo =0, (34) 

^ C ph C gr C ph c gr Cph 



with [V] as in (27), and therefore 



2c 2 2r 2 



lE osc j = — < 0, s = — > 0. (35) 

(C g r-C p h) Cph -C gr 
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In particular, there exists a two-dimensional family of phase transition waves that is parameterised by 
the speed c p h and the trivial parameter {V}. All these causality waves have no oscillations ahead of the 
interface and (SOM1) is always satisfied. The validity of (SOM2), however, depends on sgnc gr = sgnQ, 
i.e., on whether the wave is of type-I or type-II. This is illustrated by the two cartoons from Figure 4. 
We therefore find again that (SOM2) has different implications for the two wave types. This is not 
surprising since a condition on the sign of c p h — c gr does not imply any constraint for sgnc gr = sgnQ, 
or vice versa. In other words, (SOM2) does not imply the causality principle, or vice versa. For a 
type-I wave the causality principle even implies that energy is radiated towards the interface from 
behind, see the left panel in Figure 4. 

We recall that for standing sources as discussed in §3, the causality principle in fact implies (SOM2). 
It is therefore tempting to adopt the criterion (SOM2) for moving phase interfaces and reformulate 
this condition as "the energy flux has to point away from the interface, with reference to an observer 
travelling with the interface" . Let us follow this argument for causality waves as shown in Figure 4. 
Changing to the co-moving frame via j = j — c p hi and £ = £ — Cphj, the partial energy balances (9) 
transform into 

d T E osc + <9|(Q - Cph-Bosc) = +3, d T E non + dg(PV - c ph -Bnon) = -S. 

In view of Q = c gr E osc we now conclude that the causality principle implies that the relative radiation 
flux Q = Q — Cph-Eosc has the same sign as £, and is hence indeed pointing away from the interface on 
both sides. The flaw with this argument, however, is that the passage to the co-moving frame is not 
a Galilei transformation since both j and £ are Lagrangian space coordinates (Galilei transformations 
are given by Uj i— > Uj + ct). The obcrvation that Q has a certain sign on both sides of the interface 
docs therefore not imply that the interface radiates energy towards infinity. The correct treatment is 
to work with the Lagrangian space coordinates j and £ and to characterise the radiative parts of the 
energy flux in terms of Q. This flux can, as shown in Figure 4, point towards the interface from the 
tail of the wave. 

Selection rules from Riemann problems A different approach to microscopic selection criteria 
is via initial values problems. The main idea is to characterise the physically relevant solutions as the 
limit as t — > oo of reasonable initial data. For the NN chains with double well potential, a rigorous 
mathematical analysis of initial values problems is not yet available, but numerical simulations provide 
a lot of insight into the implied selection criteria. 

To illustrate this, we now present two numerical solutions of initial value problems for the NN chain 
with interaction potential (23). Due to the Hamiltonian nature of (1), we perform the simulations 
with the Verlct scheme (see, e.g., [HLW02] for details), which is a symplectic integrator and does not 
add numerical viscosity to the problem. In order to enforce the formation of macroscopic waves, we 
start with Riemann initial data, that means we set 

rj(0)=r±oo, v j (0)=v± oo for j ^ 0, 

with prescribed asymptotic data r±oo, v± OQ . Moreover, we choose a 'number of particles' N 3> 
1 and introduce the scaling parameter e = 1/N. We then integrate the lattice equation on the 
domain \j\ < N/2 over the time interval < ifi n N with < ia n < 1/2. To this end we close the 
discrete equations (2) by imposing Dirichlet boundary data for r and v. Notice that the choice of the 
computational domain and time is in accordance with both the hyperbolic scaling and Co = 1, this 
means no macroscopic wave can hit the boundary of the computational domain. 

Figures 5 and 6 depict the numerical results at time in n for two Riemann problems with sgn (r-oo) ^ 
sgn (r+oo). The left panel contains snapshots of the atomic distances against the scaled particle index 
j = sj, and the other panels show the spatial profiles of macroscopic strain R and heat flux Q, which 
are computed by mesoscopic averaging; see [DH08] for details. 
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Figure 5: Numerical solutions to the Riemann problem with = —4, r +oc = +1, v_oo = +1, 

v +OCl = — 1 and N = 4000, ts n = 0.4. The microscopic dynamics selects a type-I wave satisfying the 
causality principle. 
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Figure 6: Numerical solutions the Riemann problem with r_oo = —3/2, r +oc = +1, u_oo = +1, 
v +oa = — 1 and N = 4000, tf\ n = 0.4. The microscopic dynamics selects a type-II wave satisfying the 
causality principle. 



In both simulations we observe four elementary waves connecting five constant states. Starting from 
the left side we find: Left initial state, first wave, first non-oscillatory intermediate state, second wave, 
oscillatory intermediate state, third wave, second non-oscillatory intermediate state, fourth wave, and 
the right initial state. The first and the fourth wave, which both connect two non-oscillatory states, are 
contact discontinuities and correspond to jumps in the linear wave equation. The second wave, which 
connects to the oscillatory intermediate state, is still a harmonic phenomenon since the oscillations 
are confined to one of the quadratic wells of We also notice that the macroscopic strain does not 
change within the second wave. 

The third wave, however, describes a phase transition as it connects negative strain to positive 
strain. A closer look to the oscillatory intermediate state reveals that it is generated by a single 
harmonic travelling wave which has a certain phase velocity c p h- The third wave propagates with the 
same speed c p h, and hence it corresponds in fact to a heteroclinic travelling wave in the lattice. We 
finally mention that for e > we still have harmonic fluctuations within the constant intermediate 
state and near the contact discontinuities, but these disappear in the thermodynamic limit e — >• 0. 

Figures 5 and 6 reveal that the radiation flux Q within the oscillatory intermediate state has a 
different sign in both simulations. More precisely, we see a type-I wave in Figure 5 but a type-II wave 
in Figure 6. Recall that both waves are now 'chosen' by the microscopic dynamics (more precisely, 
by the numerical scheme). In other words, the macroscopic structure of the Riemann solution (the 
parameters of the waves and the intermediate states) determine — implicitly, but nevertheless uniquely 
- a microscopic selection rule for phase transition waves. We are not able to derive this selection 
rule rigorously from the lattice dynamics although this should, in principle, be possible. However, 
computing the speed of the interface as well as thermodynamic fields on both sides of the interface, we 
find that these meet very well the predictions from the causality principle described above, see (34) 
and (35). It would be highly desirable to understand this in greater detail and to find a reformulation 
of the causality principle that can also be applied to generic double well potentials, where Fourier 
methods can no longer be used. First steps into this direction arc already done in this paper, because 
our thermodynamic approach to radiation conditions can be generalised to more complicate nonlinear 
problems. However, we emphasise that many questions remain open. 
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Discussion The discussion in this section shows that the case of a moving interface is different from 
that of a standing source; for the latter, (SOM1) and (SOM2) are equivalent. For moving interfaces 
in phase transition waves we find that (SOM1), which is formulated in terms of sources and sinks, is 
equivalent to the entropy condition (32). The flux condition (SOM2), however, implies (SOM1) for 
type-II waves but both rules contradict each other for type-I waves. 

For moving interfaces, one has therefore to distinguish between arguments that rely on the energy 
transport in terms of fluxes, as (SOM2), and arguments based on energy productions, as (SOM1). 
Since (SOM1) is equivalent to the entropy inequality, we propose to apply this selection criterion for 
moving inhomogeneities. We also propose that flux based criteria such as (SOM2) should not be 
naively applied to problems with moving inhomogeneities. 

These suggestions are in line with the causality principle used in [SlcOl, Slc02, TV05] and numerical 
simulations of Riemann-type problems for NN chains. 

A Appendix: Macroscopic conservation laws for NN chains 

We establish the thermodynamic limit for the forced NN chain (1) provided that the forcing ( sat- 
isfies (15). Our goal is to show that the macroscopic balance laws for mass, momentum and energy 
(see (5) and (16)) can be derived rigorously as follows: 

1. The hyperbolic scaling transforms each bounded solution into a family of oscillatory functions 
which depend on the macroscopic time r and macroscopic Lagrangian space coordinate £. 

2. This family of functions is compact in the sense of Young measures, and hence we can extract con- 
vergent subsequences. Along such a subsequence, the limit measure encodes the local distribution 
functions of the oscillatory data and hence the local mean values of atomic observables. These 
local mean values provide the thermodynamic fields and are, by construction, non-oscillatory 
functions in r and £. 

3. The dynamics of NN chains implies that the thermodynamic fields of each Young measure limit 
satisfy the macroscopic conservation laws of mass, momentum and energy (in a distributional 



We now collect the mathematical tools for each of these steps. We start with some basic facts about 
Young measures and refer the reader to [Bal89, Rou97, Val94, Tay97] for more details. 

Let f2 be a domain in R fc and K be some convex and closed set in l m . A Young measure \i G 
Y(fl;K) is a f2-family of probability measures on K, that means a measurable map /i: y G CI — > 
/x(y,dQ) G Prob(K). Notice that each function Q: d — > K defines a trivial Young measure with 
fi(y,dQ) = 5Qi y \(dQ), where SQ/ y )(dQ) abbreviates the delta distribution in Q(y). 

Theorem 1 (Fundamental Theorem on Young Measures) . Each family 



is compact in the space of Young-measures YM(Cl;K). This means there exists a sequence (£„) ti( =n 
with e„ — > along with a limit measure /i G YM(Cl; K) such that 



sense) . 



(Qe) 



0<E<1 



G L°°(n;K) 



■» (-0) weakly* in L°°(fl) 



(36) 



for all observables "0 G C{K), where 




B 



gives the local mean value of tp in y G f2. 
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Proof. Sec, for instance, [Tay97], Proposition 11.3 in Section 13.11. □ 
The convergence (36) is equivalent to 

(MlO#y)dy= lim U{Q en {y))4>{y)^V (37) 

n— J-oc J 

12 n 

for all test functions <j> <G C£°(f2). Moreover, the subsequence converges strongly to some limit function 
Qo in L°°(r2;ii') if and only if the limit measure is trivial, n(y,dQ) = SQ ( y \(dQ). 

We now suppose that we are given a bounded solution to (1). As in §2, we regard the atomic 
distances rj = Mj+i — Uj and velocities Vj = iij as the basic variables, i.e., we consider 

Qj(t) = (r j (t),v j (t)), jeZ, f>0. (38) 

For a given scaling parameter < e < 1 we introduce r and £ by (4), so the macroscopic Lagrangian 
space-time coordinate is given by 

n = {(r,0 : r>0, £ G R}. 

Moreover, we identify (38) with piccewise constant functions on f2, 

Q e {et, ej + n) = Q,-(t) for every i > 0, j e Z, |t?| < §. (39) 

By assumption, we have (Qe)o< e <i C L°°(i7;if) for some ball if C K 2 , and Theorem 1 provides 
at least one subsequence that converges to some limit measure /x £ YM(Q; K). Moreover, for each 
atomistic observable ip we can compute the corresponding thermodynamic field via 

(■0(r,w))(r,£) = / i>(r,v)fj,(T,£, drdv). 
Jm. 2 

We are now able to state and prove the main result on the thermodynamic limit of forced NN chains. It 
is a direct consequence of the discrete conservation laws derived from (1), the notion of Young- measure 
convergence, and the properties of distributional derivatives. 

Theorem 2 (Macroscopic conservation laws for NN chains). The thermodynamic fields of each Young 
measure limit fi of (1) satisfy the conservation laws (5) within the following domains in the sense of 
distributions: The laws for mass (5)i and momentum (5)2 hold for tt. The conservation of energy (5)3 
holds for if £ = 0, and otherwise for Q := Q \ {£ = 0}. 

Proof. Within this proof we write Q e (t, £) = (r E (r, £), v e (t, £)) . The equation of motion (1), combined 
with the scaling rules (4) and (39), implies the following discrete conservation laws 

8 T r e - \7 +£ v e = 0, (40) 
d T v £ - V_ e *'(r e ) = C(r/e) X e(0, (41) 
a T (i« 2 + $(r e ) - eV_ e $(r £ )) - V_ 6 (« e $'(r £ )) = Cfr/eMr/e)^*) (42) 

for all r > and almost all (gl, where the discrete differential operators V±e and the scaled cut off 
function Xe ar e given by 

(v = ±fM^l±IM) , xm = wen® . 

e e 

We now multiply (40) with a test function (j> £ C£°(Q) and integrate with respect to both r and £. 
Using integration by parts and expansions with respect to e we then find 



/ r e d T 4> — v e d^(j)dTd^ = 0(e), 
Jn 
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so the limit e — >• provides (5)i in the sense of distributions, see (37). Similarly, and using that (15) 
implies 



we derive (5)2 from (41). Finally, the assertions about the energy conservation follow from (42), where 



Since the energy is conserved in we can balance the energy in the whole domain f2 via (16). 
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